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Abstract 

We present a description of symmetric nuclear matter within the framework of Landau Fermi liq- 
uid theory. The low momentum nucleon-nucleon interaction V\ ov ,^k is used to calculate the effective 
interaction between quasiparticles on the Fermi surface, from which we extract the quasiparticle 
effective mass, the nuclear compression modulus, the symmetry energy, and the anomalous orbital 
gyromagnetic ratio. The exchange of density, spin, and isospin collective excitations is included 
through the Babu-Brown induced interaction, and it is found that in the absence of three-body 
forces the self-consistent solution to the Babu-Brown equations is in poor agreement with the em- 
pirical values for the nuclear observables. This is improved by lowering the nucleon and meson 
masses according to Brown-Rho scaling, essentially by including a scalar tadpole contribution to 
the meson and nucleon masses, as well as by scaling g^. We suggest that modifying the masses 
of the exchanged mesons is equivalent to introducing a short-range three-body force, and the net 
result is that the Brown-Rho double decimation 1] is accomplished all at once. 
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I. INTRODUCTION 

Landau's theory of normal Fermi liquids P, |4| describes strongly interacting many- 
body systems in terms of weakly interacting quasiparticles. Provided that the quasiparticles 
lie sufficiently close to the Fermi surface, they will be long-lived and constitute appro- 
priate degrees of freedom for the system. The central aim of the theory is to determine 
the quasiparticle interaction, either phenomenologically or microscopically, with which it 
is possible to describe the low-energy, long-wavelength excitations of the system. This, in 
turn, is sufficient for the description of many bulk equilibrium properties of the interacting 
Fermi system. The initial application of Fermi liquid theory to nuclear physics was the 
phenomenological description of finite nuclei and nuclear matter by M igdal |(| , and later 
a microscopic approach to Fermi liquid theory based on the Brueckner-Bethe-Goldstone re- 
action matrix theory was developed by Backman Q| and others j], 0] to describe nuclear 
matter. Although the latter approach was quantitatively successful, it was observed [9] that 
Brueckner-Bethe-Goldstone theory is less reliable in the vicinity of the Fermi surface due to 
the use of angle-averaged Pauli operators and the unsymmetrical treatment of particle and 
hole self energies, which leads to an unphysical energy gap at the Fermi surface. 

With the recent development of a nearly universal low-momentum nucleon-nucleon (NN) 
interaction V[ w-fc Q| derived from renormalization group methods, the application of Fermi 
liquid theory to nuclear matter has received renewed attention ll|, [12J, [l3[ . The strong short- 
distance repulsion incorporated into all high-precision NN potential models is integrated 
out in low momentum interactions, rendering them suitable for perturbation theory calcula- 
tions. Although limited Brueckner-Hartree-Fock studies [14] indicate that saturation is not 
achieved with V\ ow -k at a fixed momentum cutoff A, it has recently been shown that by 
supplementing V[ ow _fc with the leading-order chiral three-nucleon force, nuclear matter does 
saturate, thereby justifying the use of V\ ow _k in studies of nuclear matter. 

Although many properties of the interacting ground state are beyond the scope of Fermi 
liquid theory, the quasiparticle interaction is directly related to several nuclear observables, 
including the compression modulus, symmetry energy, and anomalous orbital gyromagnetic 
ratio. As originally shown by Landau, the quasiparticle interaction is obtained from a 
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certain limit of the four-point vertex function in the particle-hole channel. It is well known 
that using realistic NN interactions in the lowest order approximation to the quasiparticle 
interaction is insufficient to stabilize nuclear matter, as evidenced by a negative value of 
the compression modulus. This general phenomenon is observed in our calculations with 
V[ ow _fc as well. However, stability is achieved by treating the exchange of density, spin, 
and isospin collective excitations to all orders in perturbation theory. The inclusion of 
these virtual collective modes in the quasiparticle interaction is carried out through the 
induced interaction formalism of Babu and Brown jl^j], which was originally developed 
for the description of liquid 3 He and later applied to nuclear matter by Sjoberg [l2, Q. 
Subsequent work [toI . 2oT ] has confirmed the importance of the induced interaction in building 
up correlations around a single quasiparticle, thereby increasing the compression modulus. 

Our study is motivated in part by the work of Schwenk et al. who were able to 

predict the spin-dependent parameters of the quasiparticle interaction from the experimen- 
tally extracted spin-independent parameters. Crucial to these calculations was a novel set 
of sum rules, derived from the induced interaction formalism, based on a similar treatment 
by Bedell and Ainsworth 2l| to liquid 3 He. In this paper we present a fully self-consistent 



solution to the Babu-Brown induced interaction equations for symmetric nuclear matter 



Our iterative solution turns out to be qualitatively similar to the results of [111 ], but we find 
that at nuclear matter density the compression modulus and symmetry energy are smaller 
than the experimentally observed values while the anomalous orbital gyromagnetic ratio is 
too large, suggesting the possibility that important phenomena have been neglected. 

We propose to extend this study by including hadronic modifications associated with the 
partial restoration of chiral symmetry at nuclear matter density, as suggested in [22]. In 
this scenario, referred to as Brown-Rho scaling, the dynamically generated hadronic masses 
drop in the approach to chiral restoration, and at nuclear matter density it is expected that 
the masses of the light hadrons (other than the masses of the pseudoscalar mesons, which 
are protected by their Goldstone nature) decrease by approximately 20%. The success of 
one-boson-exchange and chiral EFT potentials in describing the nucleon-nucleon interaction 
suggests that a modification of meson masses in medium ought to have verifiable conse- 
quences in low energy nuclear physics. Although there is much current theoretical and 
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experimental effort devoted to the program of assessing these medium modifications, the 



consequences for low-energy nuclear physics 
Applying the mass scaling suggested in 



rave yet to be fully explored. 



221 to our calculations of nuclear matter, we 



obtain a set of Fermi liquid coefficients in better agreement with both experiment and 
the nontrivial sum rules derived in [llj]. Explicit three-body forces, though essential for a 
complete description of nuclear matter, have been neglected in this study. However, we argue 
that modifying the vector meson masses is equivalent to including a specific short-ranged 
three-body force. We conclude with a discussion of the consequences of Brown-Rho scaling 
on the tensor force, which is diminished by the increasing strength of tu-meson exchange. 



II. FERMI LIQUID THEORY 

In this section we present a short description of Fermi liquid theory and its application 
to nuclear physics with emphasis on the microscopic foundation of the theory. The main 
assumption underlying Landau's description of many-body Fermi systems is that there is 
a one-to-one correspondence between states of the ideal system and states of the interact- 
ing system. As one gradually turns on the interaction, the noninteracting particles become 
"dressed" through interactions with the many-body medium and evolve into weakly interact- 
ing quasiparticles. The interacting system is in many ways similar to an ideal system in that 
the classification of energy states remains unchanged and there is a well-defined Fermi sur- 
face, but the quasiparticles acquire an effective mass m* and finite lifetimes r ~ (k — kp) -2 . 
The energy of the interacting system is a complicated functional of the quasiparticle distri- 
bution function, and in general the exact dependence is inaccessible. But one can extract 
important information about bulk properties of the system by considering small changes in 
the distribution function. Expanding to second order, one finds 

5E = ^eg^kx) + ^ E /(ki,k 2 )Mki)Mk 2 ) + 0(6n 3 ). (1) 

ki ki,k 2 

In this equation f2 is the volume of the system, ejj?' is the energy added to the system 
by introducing a single quasiparticle with momentum ki (note that for |ki| = k\ = kp, 
ej^ is just the chemical potential), and /(ki,k2) describes the interaction between two 
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quasiparticles. 

Since the quasiparticle interaction /(ki,k2) is the fundamental quantity of interest in 
Fermi liquid theory, we will carefully discuss its properties and its relationship to nuclear 
observables. Assuming the interaction to be purely exchange, it can be written as 

/(k x , k 2 ) = -rUi^k!, k 2 ) + F'(k 1; k 2 )n • r 2 + G(kx, k 2 )at ■ a 2 + G'(k u k 2 )n • r 2 a x ■ a 2 ], (2) 

iV 

where we have factored out the density of states per unit volume at the Fermi surface, 
No = ^%~t , which leaves dimensionless Fermi liquid parameters denoted by F,G,F',G'. 
The spin-orbit interaction is neglected because it vanishes in the long wavelength limit in 
which we will be interested. Also, we have not included tensor operators (which would 
greatly complicate our calculation) because the tensor force contributes almost completely 
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in second order, as shown in the original paper by Kuo and Brown 
central interaction in the 3 Si state. In the tensor Fermi liquid parameters for symmetric 
nuclear matter were calculated from Vi ow -fc in which the dominant second-order contributions 
from one-pion exchange were included. Since quasiparticles are well-defined only near the 
Fermi surface, we assume that kx — kp — k 2 . In this case the dimensionless Fermi liquid 
parameters F, F', G, G' depend on only the angle between k x and k 2 , which we call 9. Then 
it is convenient to perform a Legendre polynomial expansion as follows 

F(k,k') = £)F,P,(cos 9), G(k,k') = £G,P,(cos 9), etc. (3) 
i i 

The Fermi liquid parameters F[,Gi, . . . decrease rapidly for larger I, and so there are only a 
small number of parameters that can either be fit to experiment or calculated microscopically. 

In the original application of the theory to liquid 3 He and nuclear systems, the quasipar- 
ticle interaction was obtained phenomenologically by fitting the dimensionless Fermi liquid 
parameters to relevant data. For nuclear matter several important relationships exist be- 
tween nuclear observables and the Fermi liquid parameters. Galilean invariance can be used 
to connect the Landau parameter Fx to the quasiparticle effective mass 

m* 

— = l + Fx/3. (4) 

m 

Adding a small number of neutrons and removing the same number of protons from the 
system will increase and decrease, respectively, the density of protons and neutrons in the 
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system (and therefore the Fermi energies of the two species). The change in the energy, 
described by the symmetry energy /3, can be related j(| to the Landau parameter Fq 

h 2 k 2 

= + K)- (5) 

In a similar way, the equal increase or decrease of the proton and neutron densities leads 
to a relationship between the scalar-isoscalar Landau parameter Fq and the compression 
modulus /C 

JC = ^(1 + F ). (6) 

n 

Finally, it can be shown p that an odd nucleon added just above the Fermi sea induces a 
polarization of the medium leading to an anomalous contribution to the orbital gyromagnetic 
ratio of the form 

9? = [1 - Sgi]fi N 

9i = [$gi}nN, (7) 

where Sgi is given by 

Clearly there are certain values of the Landau parameters that are physically unreason- 
able. For instance, if F\ < —3 or F < —1, the effective mass or compression modulus would 
be negative. Quite generally it can be shown [2jJ that the Landau parameters must satisfy 
stability conditions 

X t >-(21 + 1), (9) 

where X represents F,G,F',G'. 

A rigorous foundation for the assum pti ons underlying Landau's theory can be obtained 
through formal many-body techniques j^Hl- ^ ^ s n °t our g° a l to reproduce the original 
arguments jj], but rather to give a clear motivation for the diagrammatic expansion leading 
to the quasiparticle interaction. Starting from the usual definition of the four-point Green's 
function in momentum space 

G a p, rS (ki, h; k 3 , h) = (2vr) 8 ^ 4 )(A; 1 + k 2 - k 3 - k i )[G vr (k 1 )G ps {k 2 )S W {k 1 - h) (10) 
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h + K 



fci + K 



fci + K 



FIG. 1: The Bethe-Salpeter equation for the fully irreducible vertex function V in terms of the ph 
irreducible vertex function T. 

where G{k\) is the Fourier transform of G(xt, x't') and ki, . . . , k± represent four-vectors (e.g. 
k\ = (kx, u)i)), it can be shown that the quasiparticle interaction is related to a certain limit 
of the four-point vertex function T a p jg(kii ^ 2 ; k^, £4). From energy-momentum conservation 
(k\ + ki = ks + ki) we can write k% — k\ = K = ki — k± and therefore define r(fci, k 2 \ K) = 
r(fci, k 2 \ ks, k^. The important point is that since we are considering only low-energy long- 
wavelength excitations, the particle-hole energy-momentum K should be small. We can 
write a Bethe-Salpeter equation for the fully reducible vertex function V in terms of the ph 
irreducible vertex function T in the direct channel with momentum transfer K: 

T af s^ s (ki } k 2 ;K) = r a p >jS (ki, k 2 ; K) 



shown diagrammatically in Fig. ^ The product of propagators may have singularities in the 
limit that K — > 0, in which case the poles can be replaced by 5-functions inside the integral: 



where z is the renormalization at the quasiparticle pole and 0(q) accounts for the multipair 
background. The limit K = (u, K) — > depends on the relative ordering of the two limits 
K — > and ui — > 0. Defining 




I r ae , 7r) (fc 1 , q; K)G(q)G(q + K)T^ eS (q, k 2 ; K) (11) 



G(q)G(q + K) 



2ivr^ 2 q • K 
u) — vp(\ ■ K 



6(e - fi)6(q - k F ) + 0(q), 



(12) 



T"{k x M) 

r K (h,k 2 ) 



lim lim iYfci , k 2 \ K) and 



lim lim T(ki, k 2 ; K), 



(13) 
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from eq. (fT2j) we see that the product of propagators is regular for T^. Thus, to calculate 
T w we must first calculate the ph irreducible diagrams belonging to T and then iterate via 
the Bethe-Salpeter equation with the intermediate multipair background <fi. The ^-function 
singularities in T K can be used to perform the integrals over q and |q|, and through algebraic 
manipulation it is possible to combine T w and T K into a single integral equation 

~~^ N ° z2 S / ^ r -,7#i> *&)■ ( 14 ) 

Physically, T u represents the exchange of virtual excitations between quasiparticles, and T K 
represents the forward scattering of quasiparticles at the Fermi surface. By relating these 
vertex functions to the equations describing zero sound, Landau 4] was able to make the 
identifications 

f(k h k 2 ) = z^iktM) and 

a(h,k 2 ) = z 2 T k (k x ,k 2 ), (15) 

where f(kx,k 2 ) is just the quasiparticle interaction introduced earlier and a(ki,k 2 ) is the 
physical scattering amplitude. 

III. INDUCED INTERACTION 

In principle one could exactly calculate the quasiparticle interaction by summing up all ph 
irreducible diagrams contributing to the ph vertex function in the limit k/u — > 0. Since this 
is not practicable in general, one must limit the calculation to a certain subset of diagrams. 
We could proceed by calculating the relevant diagrams order by order, but this would miss 
an essential point, which we now elaborate. From eqs. (JT4~j) and ()15j) . we see that the physical 
scattering amplitude a(ki, k 2 ) iterates the quasiparticle interaction to all orders through an 
integral equation shown schematically in Fig. El If only a finite set of diagrams are included 
in the quasiparticle interaction, then the scattering amplitude will not be antisymmetric. 
For instance, if we include only the bare particle-hole antisymmetrized vertex shown in Fig. 
Eta), then diagram (b) will be contained in the equation for the scattering amplitude but 
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FIG. 2: The diagrammatic relationship between the physical scattering amplitude a and the quasi- 
particle interaction /. 



(a) (b) (c) 

FIG. 3: Diagrams contributing to the quasiparticle interaction / and the scattering amplitude a. 
Diagrams (a) and (c) contribute to /, whereas all three contribute to a. 

its exchange diagram, labeled (c), will not. Quantitatively, the fact that the scattering 
amplitude is antisymmetric requires that it vanish in singlet-odd and triplet-odd states as 
the Landau angle 6 approaches 0. This leads to two constraints J, 27 ] on the Fermi liquid 
parameters in the form of sum rules: 

^ (l + fJ(21 + 1) + 3 1 + G\J{21 + 1)) = ° (16) 

Y(- - + - + - )=0 (17) 

Y \3 1 + ^/(2/ + 1) 1 + F//(2Z + 1)^ 1 + GV(2/+1)7 1 ' 

Clearly, the sum rules must be satisfied for the "correct" set of Fermi liquid parameters 
describing nuclear matter. To account for this infinite set of exchange diagrams, Babu and 
Brown |lj| proposed separating the quasiparticle interaction into a driving term and an 
induced term: 

f(k,k') = f d (k,V)+f i (k,k% (18) 
9 



k 2 + K 



k 2 + K 



k 2 + K k 2 




fa + K k 



1 fa + A 




+ 



fa fa + A' 




+ 



FIG. 4: The diagrammatic form of the induced interaction. In the limit that k\ = k 2 it can be 
shown that the external lines exactly couple to particle-hole excitations through the / function. 

where the induced interaction is denned to contain those diagrams that would be the ex- 
change terms necessary to preserve the antisymmetry of a(ki, k 2 ). Then the induced inter- 
action is given by a diagrammatic expansion shown in Fig. 0] Physically, the induced inter- 
action represents that part of the quasiparticle interaction that results from the exchange of 
virtual collective modes, which can be classified as density, spin, or isospin excitations. In 

ri 

the limit that ki — > h 2 it can be rigorously proved [16| that the coupling of quasiparticles to 
these collective excitations is precisely through the quasiparticle interaction itself, thereby 
justifying the diagrammatic expression in Fig. 0] 

The relationship between the induced interaction and the full quasiparticle interaction 
was derived by Babu and Brown for liquid 3 He and applied to nuclear matter by Sjoberg 



17| . To lowest order in the Fermi liquid parameters, the induced interaction is given by 



AG' 



+ 



3F 



i2 



+ 



3Grr 



+ 



9GV 2 



1 + F a 1 + F 'a 1 + G a 1 + G 'a 



+ 



3F 



i2 



Gr 



3Gq 



i2 



1 + F a 1 + F 'a Q 1 + G a 1 + G 'a 



+ 



3GV 



3G*o 



i2 



1 + F a 1 + F 'a 1 + G a 1 + Go'^o 



-^0 



+ 



G ' 2 



where 



1 + F a 1 + F 'a 1 + G a 1 + G 'a 

kp - q/2 



«0 



a 



Of, 



/ \ 1 1 / q 



q J k F + q/2 



(19) 



(20) 
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is the Lindhard function, which is related to the density-density correlation function xpp by 

-oto(q,u) 



(21) 



1 + F ao(q,c.;) 

and q = ki — k2. The interpretation of equation ()19|) is as follows. The Landau parameters 
in the numerator describe the coupling of quasiparticles to particular collective modes. For 
instance, the F represents the coupling to density excitations, Go the coupling to spin 
excitations, etc., and the denominators enter from the summation of bubbles to all orders. 
Including the I = 1 Fermi liquid parameters, the induced interaction is given by 

,2 



4F 



+ 



Fo 2 a 
1 + F a 

3Fq\q 
1 + F 'a 

where «i defined by 
ai(q,0) 



+ 1- 



q 



4A* 



+ 1- 



4^ 



Fi 2 ^ 3G 2 a _ . ( 
1 + Fxai 1 + G a 

+ z + 1 



1 + Fi'ai 1 + G 'a 



q 2 \ 


3Gi 2 ai 


4k* J 


1 + Giati 


q 2 N 


1 9G\ 2 a\ 


4& FJ 


i 1 + Gi'a 



(22) 



3 fcp 
8 ~ V 



+ 



3q 



In 



kp+_q/2 
k F - q/2 



(23) 



4g 32fc Fy 

is related to the current- current correlation function, and analogous expressions hold for the 
spin- and isospin-dependent parts of the induced interaction. These equations were first 
obtained in [20], carried far enough to include velocity-dependent effects in terms of an 
effective mass, in the approximation of quadratic spectrum. 

Having characterized the induced part of the quasiparticle interaction, let us now elabo- 
rate on the driving term. By definition, this component of the interaction consists of those 
diagrams that cannot be separated into two diagrams by cutting one particle line and one 
hole line. Some of the low order terms contributing to the driving term are shown in Fig. 
where the interaction vertices are assumed to be antisymmetrized. Some higher-order terms, 
such as diagram (d) in Fig. 03 are included implicitly through the quasiparticle renormaliza- 
tion z and need not be calculated explicitly, as described in detail in In order to preserve 
the Pauli principle sum rules ()16|) and ()17j) the driving term must be antisymmetrized. Thus, 
including Fig. Efd) requires that (e) also be included in order for the scattering amplitude 
to be antisymmetric. 
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(b) 



(c) 



(d) 



(e) 



FIG. 5: A selection of diagrams contributing to the driving term in the quasiparticle interaction. 
Diagrams (d) and (e) are included implicitly through the renormalization at the quasiparticle pole. 



IV. CALCULATIONS AND RESULTS 



According to the discussion in the previous section, the starting point of a microscopic 
derivation of the quasiparticle interaction is a calculation of the antisymmetrized driving 
term to some specified order in the bare potential. Nearly all previous calculations have used 
the G-matrix, since it is well known that the unrenormalized high-precision NN potentials 
are unsuitable for perturbation theory calculations due to the presence of a strong short- 
distance repulsion. The resummation of particle-particle ladder diagrams in the G-matrix 
softens the potential but introduces several undesirable features from the perspective of 
Fermi liquid theory. Most important is the unphysical gap in the single particle energy 
spectrum at the Fermi surface due to the fact that hole lines receive self-energy corrections 
but particle lines do not. In the past it was suggested that introducing a model space, 

within which particles and holes are treated symmetrically, could overcome this difficulty. 

An alternative method for taming the repulsive core is to integrate out the high momen- 
;um components of the interaction in such a way that the low energy dynamics are preserved 
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291 ]. This is accomplished by rewriting the half-on-sheil T-matrix 



T(p\p,p 2 ) = V NN (p\p) + -V 



7T 



VNN(p',q)T(q,p,p 2 

p2 _ g2 



-q 2 dq 



(24) 



with an explicit momentum cutoff A, which yields the low momentum T-matrix defined by 

" A Vi ow _ k (p',q)T low _ k (q,p,p 2 ) 



2 r 

Ti ow - k (p',p,p 2 ) = V low ^ k (p',p) +-V 

^ Jo 



-q dq. 



(25) 



p* — q* 

Enforcing the requirement that T low _ fc (p',p, p 2 ) = T(p',p,p 2 ) for p',p < A preserves the low 
energy physics encoded in the scattering phase shifts. Remarkably, under this construction 
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all high-precision NN potentials flow to a nearly universal low momentum interaction V[ ow -fc 
as the momentum cutoff A is lowered to 2.1 fm -1 . In fact, k = 2.1 fm -1 is precisely the CM 
momentum beyond which the experimental phase shift analysis has not been incorporated 
in the high-precision NN interactions. 

For an initial approximation to the driving term, we include the first-order antisym- 
metrized matrix element shown diagrammatically in Fig. HJa) as well as the higher order 
diagrams, such as (d) and (e), that are included implicitly through the renormalization 
strength at the quasiparticle pole. The quasiparticles are confined to a thin model space P 
near the Fermi surface 

P = hm J2 ( 26 ) 

kp<k<kp +8 

and the first-order contribution is given by 

(hhSTiviihh - hh)ST) = (k,esT\v\k,esT), (27) 

where k\ = k<i = k 3 = /c 4 = kp, 9 is the angle between the two momenta, and the relative 
momentum k = /ci?sin(#/2). Given the Vi ow -fc matrix elements in the basis \klSTJ), we 
project onto the central components and change from a spherical wave basis to a plane wave 
basis. Then the dimensionful driving term is given by 

(kST\V d \kST) = z 2 -^y2(2J+l)(l - (-l) l+s+T )(klSJT\V low ^ k \klSJT). (28) 
^ + 1 j, i 

Inserting the form of the quasiparticle interaction in eq. (J2J) into the left hand side of eq. 
(|2*H|) . we obtain the Fermi liquid parameters in terms of Vsrik) = (kST\V\kST). The result 
is 

lo lo lo lo 

1 r/ 3 1 3 

9 = -ttVoo ~ t^Voi + — Vio + —V n 
lb lo lo lo 

g . = ^-l^-l^ + lvu, (29) 
where the momentum dependence has been suppressed for simplicity. From eq. (jlj) it can 
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CD-Bonn 
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0.392 
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0.248 
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-1.199 


0.135 


0.350 


0.603 
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-0.506 


0.241 


0.252 


0.118 


-0.445 


0.161 


0.172 


0.225 


-0.498 


0.240 


0.259 


0.118 
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-0.201 


0.120 


0.101 


0.021 


-0.213 


0.127 


0.106 


0.020 


-0.200 


0.122 


0.101 


0.022 


3 


-0.110 


0.054 


0.051 


0.009 


-0.120 


0.060 


0.056 


0.007 


-0.111 


0.055 


0.051 


0.010 



TABLE I: The Fermi liquid parameters of the NN interaction V\ ow -k derived from the Nijmegen 
potentials and CD-Bonn potential for a cutoff of A = 2.1 fm _1 and Fermi momentum 1.36 fm _1 . 



be shown that 



m 
m 



T^IUS' (30) 

where \x = Imkp /7r 2 h 2 = -^Nq, from which we construct the dimensionless Fermi liquid 
parameters. In all of our calculations we include partial waves up to J = 6. In Tabled we 
show the Landau parameters of the driving term derived from three different low momentum 
interactions obtained from the Nijmegen I & II potentials j^] and the CD-Bonn potential 
3ll | for a momentum cutoff of A = 2.1 fm _1 and a Fermi momentum of kp = 1.36 fm _1 . 
From the available theoretical analyses of nucleon momentum distributions |32j , we take the 
quasiparticle renormalization strength to be z = 0.7 for nuclear matter. 

The induced interaction is obtained by iterating equations (j!8)l and (f2"2"|) until a self- 
consistent solution is reached. The density- density and current-current correlation functions 
in (}2"2"j) introduce a momentum dependence in the induced interaction, and the Fermi liq- 
uid parameters for the induced interaction are obtained by projecting onto the Legendre 
polynomials 



21 + 1 



Fi(6)Pi(cos 9)d(cos 9), etc. 



(31) 



-l 



For the first iteration we use the Landau parameters obtained from the bare low momentum 
interaction as an estimate for the full quasiparticle interaction in eq. (|22|) . However, since 
F does not satisfy the stability criteria Q for either the Nijmegen or CD-Bonn potentials, 
in the first iteration we replace it in both cases with an arbitrary value that does. The 
convergence of the iteration scheme is generally rapid and relatively insensitive to the set of 
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initial parameters chosen for the low-momentum Nijmegen I and CD-Bonn potentials. In 
contrast, the low momentum Nijmegen II potential exhibits poor convergence properties, 
though a solution to the coupled equations can still be found. For a completely consistent 
solution at each iteration we recalculate the driving term with the new effective mass. 

The final self-consistent result for the quasiparticle interaction is shown in Fig. El for 
the CD-Bonn potential, and the Fermi liquid parameters for the driving term, induced 
interaction, and full quasiparticle interaction are shown in Table ITT1 For comparison we list 




FIG. 6: The self-consistent solution for the full quasiparticle interaction as a function of k = 
||ki — k2 1 derived from the low momentum CD-Bonn potential. 

the Fermi liquid parameters obtained in where the spin-independent Landau parameters 
were taken from experiment and used to calculate the spin- dependent parameters with a set 
of nontrivial sum rules: 
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G^ 


Ft 


Gi 


F- 


Gj 





-0.476 


0.025 


0.221 


0.784 


-1.276 


0.144 


0.373 


0.642 


0.801 


-0.119 


-0.152 


0.142 


1 


-0.335 


0.263 


0.273 


0.171 


-0.530 


0.256 


0.275 


0.125 


0.195 


0.007 


-0.002 


0.048 


2 


-0.238 


0.139 


0.117 


0.020 


-0.212 


0.130 


0.107 


0.024 


-0.026 


0.009 


0.010 


-0.003 


3 


-0.101 


0.055 


0.050 


0.014 


-0.119 


0.059 


0.054 


0.011 


0.018 


-0.004 


-0.004 


0.003 



TABLE II: The self-consistent solution of the Babu-Brown equations for the low momentum CD- 
Bonn potential. The full Fermi liquid parameters are obtained by projecting the quasiparticle 
interaction in Fig. El onto the Legendre polynomials. 



F = -0.27 G = 0.15 ±0.3 F ' = 0.71 G ' = 1.0 ± 0.2 

Fi = -0.85 Gi = 0.45 ±0.3 F/ = 0.14 G x ' = 0.0 ± 0.2 

Although the experimental values for the spin-independent parameters are appreciably dif- 
ferent from the self-consistent solution we have obtained, our values for the spin- dependent 
parameters fall within the errors predicted from the sum rules. However, the main effect 
of the induced interaction is to cut down the strong attraction in the spin- independent, 
isosp in-independent part of the quasiparticle interaction. In fact, the repulsion in this chan- 
nel coming from the induced interaction is large enough for the resulting Fo to satisfy the 
stability condition in (jHJ). The effective mass, compression modulus, and symmetry energy 
are shown in Table IIHI together with the deviations 5S\ and SS^ from the sum rules (fTo|l 
and ()17|) . We list the results for the three different bare potentials with a momentum cutoff 
of A = 2.1 fm _1 . In calculating the contributions to (fho|) and (JT7j) we have included Landau 
parameters for I < 3. The compression modulus for nuclear matter is extrapolated from the 
data on g iant monopole resonances in heavy nuclei, with the expected value being 200 — 300 
MeV 0fl- The symmetry energy is determined by fitting the data on nuclear masses to 



various versions of t 
P = 25-35 MeV |3j, 



re semi-empirical mass formula |3J| , and currently the accepted value is 
. Both the compression modulus and the symmetry energy shown 
in Table IIHI are significantly smaller than the experimental values. On the other hand, the 
anomalous orbital gyromagnetic ratio, determined from giant dipole resonances in heavy 



16 





Nijmegen I 


Nijmegen II 


CD-Bonn 


m* jm 


0.887 


0.930 


0.888 


K [MeV] 


136 


102 


136 


P [MeV] 


18.1 


20.5 


17.6 




0.682 


0.452 


0.685 


SSi 


0.20 


0.16 


0.27 


5S 2 


-0.04 


-0.02 


-0.04 



TABLE III: Nuclear observables obtained from the self-consistent solution of the Babu-Brown 
equations and deviations 5 Si and 6S2 from the Pauli principle sum rules. 



nuclei, is too large compared with the experimental value of Sgf = 0.23 ± 0.03 [31 

As suggested in the introduction, we propose to remedy these discrepancies by consid- 
ering the effects of Brown-Rho scaling on hadronic masses. The proposed scaling law for 
light hadrons - other than the pseudoscalar mesons, whose masses are protected by chiral 



mvariance - is 



22, 



38] 



m 



v 



19a m N 



1-C- 



n 



(32) 

m v m a V g\ m N n 
where the subscript V denotes either the p or u vector meson, a refers to the scalar meson, 
gA is the axial vector coupling, and n/no is the ratio of the medium density to nuclear matter 
density. This scaling can be thought of as extending Walecka mean field theory, in which 
the scalar tadpole contribution to the nucleon self-energy lowers the effective mass, to the 
level of constituent quarks. Attaching a scalar tadpole on the nucleon line, as shown in Fig. 
Ufa), lowers the mass according to (|32|). and a scalar tadpole connected to the vector mesons 
gives an effective three-body force as shown in Fig.[7fb). Including the in- medium scaling of 
the axial-vector coupling, which should approach g* A — 1 at chiral restoration, the net result 
is a lowering of the in- medium m v by ~ 2/3 as much as m* N . Recent experimental results 



39 



jQ] are consistent with the scaling law (|32|) for C — 0.15 and 0.092, respectively. The 
Brown-Rho "parametric scaling" has C = 0.2. However, the dense loop term AM 4l| gives 
a shift of the p-meson pole upwards. So far no one has been able to calculate it at finite 
density. 



17 




N 



(a) 




FIG. 7: Walecka mean field contribution of the scalar tadpole to the nucleon mass (a) and its 
extrapolation to constituent quarks in vector mesons (b). 

A number of previous studies (42, were successful in describing nuclear matter 

by starting from a chiral Lagrangian with nuclear, scalar, and vector degrees of freedom in 
which the hadronic masses were scaled with density according to (J32j) . In particular, the 



compression modulus and anomalous orbital gyromagnetic ratio were found to be in excellent 
agreement with experiment, which suggests that a similar approach may prove fruitful in 
our present analysis. An alternative approach, complementary to the chiral Lagrangian 
method, is to include medium modifications directly into a one-boson-exchange potential. 
Such a calculation was carried out in J3] to study the saturation of nuclear matter. In 
their work it was suggested that the a particle should be constructed microscopically as a 
pair of correlated pions interacting largely through crossed-channel p exchange. Medium 
modifications to the a mass then arises naturally from the density-dependence of the p 
mass. The final conclusion established in j3] is that at low densities the a scales according 
to (J32j) but that toward nuclear matter density the scaling is slowed to such an extent that 
saturation can be achieved. 

We proceed along the lines of Q| and introduce medium modifications directly into a one- 
boson-exchange potential. The most refined NN potentials in this category are the Nijmegen 
I, Nijmegen II, and CD-Bonn potentials. The Nijmegen potentials include contributions from 
the exchange of p, u, <fi, a, fo, and cto mesons, as well as the pseudoscalar particles which do 
not receive medium modifications in Brown-Rho scaling. The CD-Bonn potential includes 
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two vector particles (the p and u) and two scalars (<7i and 02 )■ For both potentials we 
scale the vector meson masses by 15% and the scalar meson masses by 7%. In this way we 
roughly account for the decreased scaling of the scalar particle mass observed in 

urn- 

In the 

full many-body calculation we also scale the nucleon mass by 15% and with an additional 
a/ g* A j gA — 1/ Vl-25 at nuclear matter density It is essential to also scale the form factor 
cutoffs Af of the vector mesons in the boson-exchange potentials. 

In Table IIVI we show the effective mass, compression modulus, symmetry energy, and 
anomalous orbital gyromagnetic ratio for the Nijmegen I & II and CD-Bonn potentials 
with the in-medium modifications. We also show for comparison the results from the Ni- 
jmegen93 one-boson-exchange potential, which has only 15 free parameters and is not fine- 
tuned separately in each partial wave. We observe that the iterative solution is in better 
agreement with all nuclear observables. The anomalously large compression modulus in 
the CD-Bonn potential results almost completely from the presence of a large uj coup ling 
constant g^^/A-K = 20.0. With the same g^NN/^ anc ^ Bonn-B potential, Rapp et al. |45j 
obtained /C = 356 MeV. The compression modulus is very sensitive to this parameter, as we 
have checked that dropping this coupling by 20% cuts the compression modulus in half but 
alters the other nuclear observables by less than 5%. The naive quark model predicts a ratio 
°f 9^nn 1 9pNN = 9 between the u and p coupling constants, which is largely violated in the 
CD-Bonn potential g^ NN / 'g 2 p NN = 24 though roughly satisfied in the Nijmegen potentials 
QuNn/ 9pNN = 11' perhaps resulting in better agreement with experiment. 

Thus, by extension of the Walecka mean field on nucleons to those on constituent quarks, 
we obtain the Fermi liquid parameters for the theory that is now essentially Brown-Rho 
scaled, as shown in Table One should note that these results are only for infinite nuclear 
matter and especially the three-body term will act in many different diagrams in the finite 
systems. However, our arguments suggest that the three-body terms intrinsic to Brown-Rho 
scaling will be useful in stabilizing light nuclei. 
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V m 


^NII 


^V93 


VcDB 


m*/m 


0.721 


0.763 


0.696 


0.682 


K [MeV] 


218 


142 


190 


495 


P [MeV] 


20.4 


25.5 


23.7 


19.2 


$91 


0.246 


0.181 


0.283 


0.267 



TABLE IV: Nuclear observables obtained from the self-consistent solution to the Babu-Brown 
equations incorporating Brown- Rho scaling. Four different bare potentials - the CD-Bonn potential 
(Vcdb), Nijmegen I (Vni), Nijmegen II (Vnii), and Nijmegen93 (iV93) potentials - were used to 
construct low momentum interactions for a cutoff of A = 2.1 fm _1 . In eq. (|32|) the parameter 
C = 0.15. 



1 


Fi 


Gi 




G\ 





-0.20 ± 0.39 


0.04 ± 0.11 


0.24 ± 0.16 


0.53 ± 0.09 


1 


-0.86 ± 0.10 


0.19 ± 0.06 


0.18 ± 0.05 


0.17 ± 0.12 


2 


-0.21 ± 0.01 


0.12 ± 0.01 


0.10 ± 0.02 


0.01 ± 0.02 


3 


-0.09 ± 0.01 


0.05 ± 0.01 


0.05 ± 0.01 


0.01 ± 0.01 



TABLE V: Fermi liquid coefficients for the self-consistent solution to the Babu-Brown equations 
using Brown-Rho scaled nucleon and meson masses in the four low momentum CD-Bonn and 
Nijmegen potentials listed in Table ITVl The tabulated values display the average and spread from 
the four different potentials and not the actual uncertainties associated with the Fermi liquid 
parameters. 

V. DISCUSSION OF THE TENSOR FORCE WITH DROPPING p-MASS IN SAT- 
URATION OF NUCLEAR MATTER 

The tensor force contributes chiefly in second order perturbation theory as an effective 
central force in the 1 = 1 channel. As the density increases, some of the intermediate states 
are blocked by the Pauli principle. In the two-body system the tensor force contributes to 
the 3 Si state, but not to the 1 S'o state, and gives most of the attractive interaction difference 
between the 3 5*i and 1 S states, effectively binding the deuteron. However, the intermediate 
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FIG. 8: Reduction in the strength of the tensor force due to a scaled /9-meson mass. Contributions 
from both 7r-meson and p-meson exchange are included in both curves. We have used the Brown- 
Rho parametric scaling, so that at nuclear matter density m* = 0.8m p . 

state energies relevant for the second-order tensor force are > 225 MeV (See Fig. 69 of |3| 
which is for 40 Ca. For nuclear matter the intermediate state momenta would be higher.), 
well above the Fermi energy of nuclear matter, and most intermediate momenta are above 
the Viow-fc upper model space limit of 420 MeV/c, so the tensor force is largely integrated 
out. 

However, since the beginning of Brown- Rho scaling it has been understood that the tensor 
force is rapidly cut down with increasing density. That is because the pion mass does not 
change with density, being protected by chiral invariance, but the p-meson mass, which is 
dynamically generated, decreases by 20% (parametric scaling) in going from a density of 
n = to nuclear matter density n — no- Since the p- meson exchange contributes with 
opposite sign from that of the pion, this cuts down the tensor force substantially. In Fig. |H] 
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we show the total tensor force from ir and p exchange at zero density and nuclear matter 
density no- Since it enters in the square, this means a factor of several drop in the tensor 
contribution to the binding energy, as shown in Fig. El 

We believe that the work of ref. [3^ shows unambiguously that the mass of the cu-meson 
is ~ 14% lower at nuclear matter density than in free space. It is remarkable that nu- 
clear structure calculations have been carried out for many years without density-dependent 
masses but with results usually in quantitative agreement with experiment. In Jj Brown 
and Rho showed that in cases where the exchange of the 7r-meson is not important, such as 
in Dirac phenomenology, there is a scale invariance such that if the masses of all relevant 
mesons are changed by the same amount, the results for the physical phenomena are very 
little changed. 

Since the pion exchange gives the longest range part of the nucleon-nucleon interaction, 
it is amazing that there are not clearcut examples in nuclear spectroscopy such as level 
orderings that are altered by the p-meson exchange playing counterpoint to the 7r-meson 
exchange, as we find in this paper for nuclear saturation. The in-medium decrease in the 
p-mass increases the effect of p-exchange, which enters so as to cut down the overall tensor 
force, the p and n exchange entering with opposite sign. 

Finally, nearly forty years since the Kuo-Brown nucleon-nucleon forces were first pub- 



lished, it was shown 



that the summation of core polarization diagrams to all orders 



is well-approximated by a single bubble. However, in light of the double decimation of 
being carried out here in one step, these forces should be modified to include the medium 
dependence of the masses. Phenomenologically this can be done by introducing three-body 
terms, as we did here, but from our treatment of the second-order tensor force it is clear 
that this should be done at constituent quark level. 

VI. CONCLUSION 

We believe that by discussing the nuclear many-body problem within the context o ' 
Fermi liquid theory with the interaction Vi ow _fc following the work of Schwenk et al. [ll| 
we have a format for understanding connections between the physical properties of the 
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10 




r [fm] 



FIG. 9: Reduction of the tensor force in second order perturbation theory due to a scaled p-meson 
mass. The intermediate state energy is approximated as 225 MeV. Contributions from both tt- 
meson and /9-meson exchange are included in both curves. At nuclear matter density, no, we have 
used the parametric scaling m* = 0.8m p . 

many-body system and the nuclear potentials. We carried out an iterative solution of the 
Babu-Brown equations, which include both density-density and current-current correlation 
functions, calculating input potentials via a momentum space decimation to V\ ow -k- By 
including Brown-Rho scaling through scalar tadpoles, as suggested by Walecka theory, our 
iterative solution provides the empirical Fermi liquid quantities. Our nucleon effective mass 
is on the low side of those usually employed, as is common in Walecka mean field theory. 
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